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Compact and extended dendrimers are two important classes of dendritic polymers. The impact of the 
underlying structure of compact dendrimers on dynamical processes has been much studied, yet the relation 
between the dynamical and structural properties of extended dendrimers remains not well understood. In this 
paper, we study the trapping problem in extended dendrimers with generation-dependent segment lengths, 
which is different from that of compact dendrimers where the length of the linear segments is fixed. We first 
consider a particular case that the deep trap is located at the central node, and derive an exact formula 
for the average trapping time (ATT) defined as the average of the source-to-trap mean first passage time 
over all starting points. Then, using the obtained result we deduce a closed-form expression for the ATT 
to an arbitrary trap node, based on which we further obtain an explicit solution to the ATT corresponding 
to the trapping issue with the trap uniformly distributed in the polymer systems. We show that the trap 
location has a substantial influence on the trapping efficiency measured by the ATT, which increases with the 
shortest distance from the trap to the central node, a phenomenon similar to that for compact dendrimers. 
In contrast to this resemblance, the leading terms of ATTs for the three trapping problems differ drastically 
between extended and compact dendrimers, with the trapping processes in the extended dendrimers being 
less efficient than in compact dendrimers. 

TAGS numbers: 36.20.-r, 05.40.Fb, 05.60.Cd 


I. INTRODUCTION 

As an integral class of macromoleculesi, dendrimers 
have attracted extensive attention of scientists in the in¬ 
terdisciplinary fields of physics, chemistry, and materi¬ 
als, since they display unusual geometrical, physical, and 
chemical properties^Ti^. These chemical compounds are 
synthesized by repeating units arranged in a hierarchical 
self-similar fashion, and are characterized by basic build¬ 
ing element, branching of the end groups, as well as the 
number of generations. Among various dendritic poly¬ 
mers, compact and extended dendrimers constitute two 
important types of dendrimer families. In the compact 
family, the length of linear unit (segment) in all genera¬ 
tions is identical, while in the extended family segment 
length is varying, which decreases toward the periphery. 
The unique structure of two dendrimers makes them po¬ 
tentially promising candidates for a broad ranges of ap¬ 
plications, e.g., artificial antenna systems for light har¬ 
vesting^!^. 

Light harvesting by dendrimers can be described by 
trapping process^^— , which is a kind of random walks 
with a deep trap positioned at a given site, absorbing all 
particles visiting it. The highly desirable quantity for this 
paradigmatic dynamical process is trapping time (TT), 
also known as mean first-passage time (MFPT)i^“— , 
which is the expected time for a particle starting off from 
a source point to first reach the trap. The average trap- 
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ping time (ATT), defined as the average of trapping time 
over all starting nodes, gives insight to the trapping pro¬ 
cess, since it provides a quantitative measure of trapping 
efficiency. In addition to light harvesting, trapping is 
also closely related to many other significant dynami¬ 
cal processes in diverse complex systems, such as target 
searcbi^— , energy or exciton transport in polymer sys- 
tems2£“— . 

In consideration of its direct relevance, trapping in 
various complex systems has been an active subject of 
research in the past years, with particular attention fo¬ 
cused on determining ATT in diverse systems that dis¬ 
play different structural and other properties. Thus 
far, trapping problem has been extensively studied for 
many systems, including regular planar and cubic lat- 
tices^^— , T— fractals and their extensions'^— , Sierpin- 
ski gasket^2ii^ and Sierpinski tower—, small-world uni¬ 
form recursive trees^^— , scale-free networks^^— , com¬ 
plex weighted networks^^— , as well as directed weighted 
treelike fractals^. These studies uncovered how the be¬ 
havior of ATT for trapping is affected by different prop¬ 
erties of complex systems, such as topological structure, 
weight and direction of edges. 

Except for the aforementioned systems, concerted the¬ 
oretical efforts have also been devoted to trapping prob¬ 
lem in polymer systems (especially dendrimers), report¬ 
ing specific properties and phenomena of the trapping 
process in these macromolecules. The trapping problem 
in compact dendrimers was first addressed in Refsi2r2i^ 
where the MFPT from the peripheral node to the cen¬ 
tral node was computed both in the presence and in the 
absence of a fixed energy bias. Partly inspired by these 
works, applying the theory of finite Markov chain, the ex- 
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act expression for MFPT from an arbitrary node to the 
central node was derived in Refsi^^— , where the ATT 
was also obtained when an immobile trap is located at the 
center. The influence of trap location on ATT for trap¬ 
ping in compact dendrimers was explored in^. Moreover, 
the factors governing the trapping efficiency for compact 
dendrimers have also been much studied in detail^^— . 

In contrast to compact dendrimers, trapping problem 
in extended dendrimers has received less attention^, al¬ 
though they are an important class of nanoscale super¬ 
molecules and their specific structure is suggested to yield 
different trapping behavior from that for their compact 
counterparts. Since extended dendrimers are an impor¬ 
tant class of nanoscale supermolecules, it is of theoreti¬ 
cal and practical importance to address the influences of 
topology of extended dendrimers on trapping occurring 
on them. 

In this paper, we present a comprehensive study of 
trapping in extended dendrimers, with an aim to explore 
the impact of their geometrical structure on the trapping 
efficiency. We first focus on a special case of trapping 
problem with a single trap located at the central node. 
We derive closed-form expressions for the MFPT from 
an arbitrary node to the trap, as well as for the ATT 
to the trap over all starting points. Then, using these 
results, we determine an exact solution for the ATT to 
an arbitrary target. Finally, we attack the case of trap¬ 
ping with the trap distributed uniformly over all nodes. 
By using two different techniques, we deduce an explicit 
formula for the ATT. We show that the trap position has 
a strong effect on the efficiency of trapping occurring on 
extended dendrimers, since their leading scalings display 
distinct dependence on the system size, which increase 
with the shortest distance from the trap to the central 
node. We also demonstrate that for all the cases of trap¬ 
ping problems, the dominant behaviors for trapping in 
extended dendrimers are different from those correspond¬ 
ing to compact dendrimers. 


II. CONSTRUCTION AND PROPERTIES 

In this section, we introduce the construction of ex¬ 
tended dendrimers and study their relevant properties. 


A. Construction method 

We concentrate on a particular network for extended 
dendrimers, with the number of nearest neighbors of any 
branching site being 3. Let Dg (g > I) denote the ex¬ 
tended dendrimer having generation number g, with the 
length of segment at generation i being g — i. Then, 
Dg can be built by g iterations in the following wa y^^’^° . 
At iteration 0, the dendrimer comprises only one central 
node; at iteration I, 3 segments with length g — 1 are 
created being attached to the central node. At iteration 
i {1 < i < g), for each branching node of generation 
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i — 1 (i.e., the node generated at iteration i — 1 farther¬ 
most from the center), 2 segments (chains) containing 
g — i nodes are generated and are linked to the branching 
node. At the last iteration g, for each node created at 
generation g — 1 , we add a pair of new nodes and link 
them to it. All the nodes introduced at this stage are 
called peripheral nodes of Dg. Figure [T] illustrates the 
structure for a specific extended dendrimer D 4 . 

Let Mg denote the length of the shortest path from an 
arbitrary peripheral node to the center. Then, 

s-i . 

i=l 

According to their shortest distances to the central node, 
all nodes in Dg can be divided into Mg + \ levels, with the 
central node being at level 0 , and the peripheral nodes 
being at level Mg. 

Let Gi{g) denote the generation (iteration) at which 
all nodes at level i were created. In order to determine 
Gi{g), we introduce two quantities Li{g) and Si, where 
Li (g) represents the possible maximal value of levels that 
nodes created at generation i belong to, and Si is defined 
to be 


Sr 


1, Z = 0, 

|i + Ej, * > 1, 

+ i+ 2 

2 ■ 


Then, we have 


Lr{g) 


Mg-Sg-i-i, 0 < i < g, 
Mg, i = g. 


( 2 ) 


(3) 
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and 


2. Average path length 


G^{9) = 


0 , 

J I J S ~ Sg-j,Mg — Sg-j-l], 

9, 


i = 0, 

0 < l < Mg, 
i = Mg. 

(4) 


By construction, the number of nodes at level i is 


The average path length denotes the average of length 
for the shortest path between two nodes over all node 
pairs. Assume that each edge in Dg has a unit length. 
Then the length of the shortest path from node i to j in 
Dg, denoted by dij{ 9 ), is the minimum length of the path 
connecting the two nodes. Let dg represent the average 
path length of Dg, defined by 


N^{9) 


1 , 1 = 0 , 

3.2Gds)-i^ l<i<Mg. 


Thus, the total number of nodes in Dg is 


(5) 


T _ ^tot(g) 

® " Ng{Ng - l)/2 ’ 


( 11 ) 


where Ttot(g) is the sum of dij{g) over all pairs of nodes, 
i.e.. 


= =9-29-i-3g-2. ( 6 ) 

i^O 

And the total number of edges in Dg is 

Eg = Ng-1 = 9- - 3g - 3. (7) 


B. Structural properties 

We proceed to present some important structural prop¬ 
erties of Dg, focusing on degree distribution and average 
path length. 


1. Degree distribution 

It is obvious that, for nodes in Dg, there are only three 
possible degree values (1, 2, and 3). The degree of any 
node at levels Lo{g), Li{g), ..., Lg-i{g) is 3; the degree 
of all peripheral nodes (i.e., nodes at level Mg) is 1; and 
the degree of all other nodes is 2. Let Ai{g) represent 
the number of nodes with degree i. Then, it is easy to 
obtain 

Ai(g) = iVM,(g) = 3-2®-i, ( 8 ) 

g-2 Li{g)-1 

^ 2 ( 9 ) = 

i=l j=Li^i(g) + l 

= 3 • 29 - 3 • 29-1 -39 + 1 , (9) 

and 

3-1 

A3(g)=^iVL,(,)= 3-29-1-3, 

2=1 


Ptotig) — dij{g). (12) 

i<j 

We continue by showing the procedure of determin¬ 
ing Ptot(g), which just equals the number of edges in the 
shortest paths between all pairs of nodes in Dg. Instead 
of counting the edges in the paths, here we count the 
paths passing through a given edge, and then sum the 
results of all edges in Dg. Let {i,j) be an edge in Dg 
linking two nodes i and j, and let Eij(g) be the number 
of the shortest paths of different node pairs, which pass 
through {i,j). Let Ni^j{g) denote the number of nodes 
in Dg lying closer to node i than to node j, including i it¬ 
self. Then, total shortest-path distance can be computed 

throug h I’ld^ 

Ptot{g) — ^ ^ Eijig) — ^ ^ dVi<j{ 9 )Nj<i{g) 

{i,j)eDg {i,j)eDg 

— 'D/ ^i<ji9)[d^g ~ ^i<ji9)] J ( 13 ) 

{i,3)eDg 

where we have used the relation Nj^i (g) = Ng — Ni^j {g). 

We now employ the relation in Eq. (USD to derive 
Ptot{g)- To this end, we look upon Dg as a rooted tree 
with the central node being the root, and use Bi{g) to 
denote the number of nodes in the subtree with a node 
at level i being its root. Then, it follows that 

Biig) = ~ 1 ) + pGi(g){9) — * + 1 

= gG,{g) + 3 • - 2g - i - ^G,{g)^ 

+ 2^*(ff) “ 1- (14) 

For any edge {i,j) connecting two nodes i and j, we 
assume that node j is the father of node i in the rooted 
tree, and that node i is at level 1. Then, for the two 
adjacent nodes i and j, Ni^j{g) = Bi{g) always holds. 
Combining this fact and the results given by Eqs. (USD 
and (flTl) . the quantity Ptot{g) can be evaluated as 

Mg 

Ptotig) = N^ig)B,{g)[Ng - B^ig)] 

2=1 


respectively. 
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= i (si • 49^2 - 189 ■ 4^g - 189 ■ 2^®+^ - 9 • 2®+^/ 
8 V 

+ 45 • 29+2^2 44J.29+1^ _ 105.29+2 _ 52^3 

-120g2+ 1243 + 816). (15) 

Inserting Eq. (fT^ into Eq. dm gives 

“ 8 (81 • 229-2 _ 27 • 2ag - 45 • 29-^ + 9^2 + 153 + 6 ) ' 

(162 • 4932 - 378 • 493 - 189 • 229+2 _ 9 .29+2^3 
+45 ■ 29+332 + 441 • 29+23 - 105 • 29+3 _ 404^3 
-2403^ + 2483+ 1632). (16) 

Equation ([6|) shows that in limit of large network, 


3 - log 2 Ng , 

which, together with and Eq. (HU), yields 
162•4932 


"" 8 • 81 ■ 229-2 


(J^NgY 


(17) 


(18) 


for very large systems. Thus, extended dendrimers are 
not small-world systems, which is in contrast to compact 
dendrimers that display the small-world phenomenons^ 
with their average path length growing logarithmically 
with the system size. 

Actually, in addition to the average path length, the 
diameter of Dg, denoted by ftg, also scales with system 
size Ng as fig ^ (In NgY, which can be obtained from the 
following arguments. For a networked system, its diam¬ 
eter is defined as the maximum of the shortest distances 
between all pairs of nodes in the system. By construc¬ 
tion, it is easy to prove that fig = 2 Mg = 32 — 3 + 2 , 
which displays a similar scaling as that of the average 
path length for large systems. 


III. TRAPPING IN EXTENDED DENDRIMERS WITH A 
PERFECT TRAP 

After introducing the construction and structural 
properties of extended dendrimers, in this section we 
study the trapping problem in Dg, in order to gain a 
comprehensive understanding of trapping process on this 
dendrimer family and explore the influence of their inter¬ 
nal structure on the trapping efficiency. We concentrate 
on three cases of trapping problems. In the first case, 
we consider trapping with a single trap positioned on 
the central node; for the second case, we attack trapping 
with the immobile trap at an arbitrary node; and in the 
last case, we address trapping with the trap distributed 
uniformly. 

The trapping problem studied here is a kind of 
isotropic discrete-time random walks in Dg with a sin¬ 
gle trap. At each time step, the walker jumps from its 
current location to an arbitrary nearest neighbor with 
identical probability. Let Tij{g) denote the MFPT from 


node i to j. Thus, Tij{g) = 0 for f = j. Let Tj{g) be the 
ATT to trap node j. Then, this interesting quantity is 
given by 

= (19) 

9 i 

Next we will study Tj{g) for the three cases of trapping 
problems defined in Dg. 


A. Trapping with a deep trap at the central node 

We first investigate the case that the trap is located 
at the central node in Dg. By symmetry, all nodes at 
the same level have an identical TT to the trap. For the 
sake of brevity, we use Pi( 3 ) to represent the TT for a 
node at level i in Dg. Then, Fi{g) satisfies the following 
relations: 


fo. 


i = 0, 


Fiig) = 


liFi-iig) + 1 ] 

+§[-^*-1-1(5) + i]i 


^[Fi-i{g) + 1 ] 

+^[Fi+i{g) + 1 ], 


0 < f < Mg and i = La.(^g){g), 


0 < f < Mg and i Y LGi{g){g), 


.pMg-i(3) + l, i = Mg. 

( 20 ) 

For f = 0 and i = Mg, Eq. (1201) is obvious; while 
for 0 < i < Mg, Eq. (|20)) can be explained as follows. 
Let Pi and p 2 represent the probabilities for a walker 
starting from a node at level i and taking one time step 
to arrive at a neighboring node at level i — 1 and level 
i + I, respectively. Since the walker performs isotropic 
random walks, for i = Lg.(g)( 3 ), pi = ^ and P 2 = |; and 
Pi = P 2 = ^ otherwise. For 0 < f < Mg, the first term 
on the right-hand side describes the process that with 
probability pi the walker starting from a node at level i 
takes one step to arrive at its unique neighbor at level 
i — I, and then makes Fi-i{g) jumps to reach the trap 
for the first time; the second term accounts for the fact 
that with probability p 2 the walker first makes a jump at 
a neighbor at level i + 1 and then takes P)_|_i ( 3 ) steps to 
visit the trap for the first time. 

From Eq. (l20l) . we can derive that for 0 < i < Mg, 


Fr{g) - F,_i{g) = 2 [P)+i(3 ) - +,( 3 )] + 3 (21) 


and 


FYg) - F,_i(3) = F,+i(3) - FYg) + 2 (22) 


hold for i = LGi(g){g) and i Y FGi{g){g), respectively. 
Let Ai{g) = Fi{g) - Fi-i{g). Then, 


A{g) 


2A,+i{g) + 3, i = LG,[g){g), 

^*- 1 - 1 ( 5 )+ 2 , i Y LGi(g){.g), 


(23) 
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is true for 0 < i < Mg. Considering AMg[g) = FMg{g) — 
FMg-i{g) = 1, Eq. (1^ can be solved to yield 


Mg) = 3.2-G^(9)+s+i-G,(g)2 + ( 25 + 3 )G,(g)- 2 i- 4 g- 3 . 

(24) 

Applying the obtained intermediate quantity Ai{g), 
the TT from a node at level i (0 < i < Mg) to the 
central node can be calculated by 


F^ig) = Fo{g) + [fM) - Fj-iig) 

i 

= Fo{g) +'^Aj{g) 
i=i 

= 2 ^)F'i{g)^ — ^(^5^ + 8g — l)Gi{g)^ 

+ i(4g2 + 4g- 3)G,{g) + (12 • 2®+^ - 3 • 2S+^g) ■ 

^-{Gi{g)+2)Q^^g^ + (6 - 

-Gi{g)^i + {2g + 3)G,(g)i + 3 • 

- 4(g + l)i + 3(23+^g - 25+2). 

(25) 

Then, according to Eqs. dm) and (l25|) . the closed-form 
expression for ATT to the trap at the central node in Dg^ 
denoted by Tc{g), can be obtained as 


Mg 

E NM)FM) 

1 


; ^27 • 4S • g - 27 ■ 2^®+^ 


9-29-1-35-3 
-9 • 29^3 - 9 • 29-2 . ^2 45 .23-2 . ^ 099.23-i 




(26) 


For a large system, i.e., g —>■ oo, Tc{g) has the following 
dominant term; 


27 • 49 ■ u 

rc(s) = -g: 2 P#~ 2 »s, 


(27) 


from which we can obtain the dependence relation of 
Tc{g) on the network size Ng as 


Tc{g)^ NglnNg, (28) 


a scaling different from that for trapping in compact den- 
drimers, where the ATT to the center node is a linear 
function of the network size^i^. 


B. Trapping with the trap at an arbitrary node 

In Sec. IIII A[ we obtained the ATT to the central node 
in Dg. Here, we use this result to further study another 


trapping problem with the trap located at an arbitrary 
node. For convenience of the following description, let 
ri represent a node at level i. Notice that all nodes at 
the same level are equivalent to one another, in the sense 
that the ATT is the same, if any of them is considered as 
a trap. Let T™™( 5 ) denote the sum of the MFPT from 
a starting point to a target node at level i {0 < i < g) 
in Dg, where the sum runs over all starting nodes in Dg. 
That is. 


TM(9)= E (29) 

jeDg 


Then, the ATT to an arbitrary node at level i in Dg is 

TrM) = j;^TMig)- (30) 

Thus, to find Tr-{g), we can alternatively evaluate the 
quantity T®“™( 5 ). 

To determine TM^(g)y consider Dg as a rooted tree 
with its root being the central node. Then, we have the 
following relation 


TM(g) 

= + 


{Ng - 1) + 1 


2-1 

Froriig) F E! Frin (g) 
i=i 


Gi{g)-1 

- E 


^T^rig - j) + l{Ng-j - 


_|__j.sum(^ - Gi{g)) + -^{^g-Gi{g) “ 1 ) ' 

LGi{g){a) 

[FLa,ig){g){9) - F^{g)] + E E'(£') “ FM)], 

j=i+l 


(31) 


where T^-n {j < i) is the MFPT from a node at level j 
to one of its descendant node at level i, Fj{g) (or Fi{g)) 
denotes the MFPT from a node at level j (or i) to the 
central node, and Fj{g) — Fi{g) {j > i) is thus equal to 
the MFPT from a node at level j to its ancestor node at 
level i. 

Equation (I5T]) can be elaborated as follows. The first 
two terms describe the process that a walker starting off 
from a node, which and the trap node have the central 
node as the lowest common ancestor, should first visit 
the central node, and then jumps (g) more steps to 
hit the trap for the first time. Here, the lowest com¬ 
mon ancestor for two nodes i and j is the node with the 
possible biggest level value in the rooted tree but having 
both i and j as its offspring nodes. The third and forth 
terms explain the sum of MFPTs to the trap, running 
over those starting nodes, which and the trap have the 
lowest common ancestor at level j(0 < j < i). The re¬ 
maining terms account for the sum of MFPTs from all 
descendants of the trap to the trap itself. 
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We now deduce the two quantities T^™(g) and Tr^r^ (ff) 
with i < j. First, making use of Eq. (I^ . we can easily 
determine T^^™{g) given by 

T^ria) = Tcig) • (Ng -1 ) 

= 1 (27 • 43 +I 5 - 27 • 2^3+3 - 9 • -9-2Sg^ 

+45 • 29g + 699 • 23+1 _ _ 294^2 _ 326 ^ 

-1176). 

(32) 

We continue to determine Tr^rjig) {i < j), which can 
be expressed as 

1-1 

Tr,rj ( 5 ) = X! (s) • (33) 

k—i 

Before determining Tr^rj (g), we first calculate the MFPT, 
Tr^n+iig), from a node at level i (0 < i < Pig) to its 
neighboring node at level i + 1. For i = 0, we have 

Troriig) = 3 + 3 [l + -F"! (ff) + Trgn (ff)] , (34) 

which, using Eq. (l25)l . is solved to yield 

T,„,,(5) =3-23+1-4g-3. ( 35 ) 

For 0 < i < Pig, the following relations can be estab¬ 
lished: 

TriTi^l (g) 

3 + 3 [l + Tri-iri [q) + (g)] 

+-F'i(fl — *) + i = La^{g)ig), 

, '5 + I (g) + (s')]’ * ^Giig)ig)- 

(36) 

Considering the initial condition in Eq. (l35l) . Eq. (l36l) is 

I 


solved to give 


'9-23-45-3-23“iid9) 

-2Gi{g)-3, i = LGi(g){g), 


Tr 


i^i + l 


9 ■ 2® - 3 - 2p 
_3.29-G.(9)+i_2gGi(g) 
PGi{g)^ — 3Gi{g) + 2t — 1, 


(37) 


* 7^ LCi{g){g)- 


Substituting Eq. (| 37 l) into Eq. (l 33 l) leads to 
1-1 

TriTj (s) = Trurk+l (s) 
k—i 

= (3 • 23 -Gd 9 ) _ g 2 _ 2^ ^ 

- 3 g ■ 23-G.(9)+i + 9.29-Gds) +5-0 Gi{g) 

+ 02 - 3 • 23 -Gi( 9 ) + 25-0 Gg{g)^ + (^ - 5" 

+35 • 23-Gds)+i _ 9.2S-Gi(9) _ g + 0 Gj(5) 

+ [(25 + 3)G,(5) + 3 • 23-Gd9)+i _ 9 .29 + 25 

—Gi{g)^ + 2 ]z + 0+2) Gi{g)^ + [ — { 2 g + 3 )Gj{g) 
-3 • 23-^1 (3)-+ + 9 .29 _ 2 g + Gj{g)^ - 2 ]j 

+3 • 23-Gd9)-H2 _ 3^ . 2S-G,(g)-Hl _ ^ 0 Q^^g)^ 

-3 • 23“'^i( 3)-I-2 _)_ 3^ . 29-G3(g)-l-l _ j2 _|_ j2 

-\G,ig)^ + \Ggig)\ 

( 38 ) 


Finally, plugging Eqs. (1^ . (15^ . and (1551) into 
Eq. (1511) , then inserting the obtained result into Eq. (1501) , 
we analytically obtain a rigorous expression for the ATT 
to any trap node at level i of the extended dendrimer 

Dg: 


TrM = 


T?rig) 

Ng-l 


1 

9 • 29-1 -Sg-S ■ 

[( - 9 ■ 2® - 6Gi(5)2 + 6Gi(5) - 14)5® + (9 • 2^Gi{g)^ - 9 • 2^Gi{g) - 9Gi(5)2®"‘^‘(»^+2 - 9 • + 135 • 23"^ 

+12iGi(5) - 61 + 6Gi(5)^ - 17Gi(5)2 + llGi(5) - ^)g^ + ( - 9 • 23+iiGi(5) + - 9 • 23+ii 

-9 • 2‘>Gi(g)^ + 9 • 23+iGi(5)2 + 9Gi{g)^20-^'^^'>+^ - 9 • 2‘>Gi{g) + 3Gi{g)20-^'^o'>+^ + 27Gi{g)2'^o-^'^o'>+^ 

+21 • 2s-iih9)+i 27 ■ 229-Gd9)+i _ ;^23 • 2®"^ - 27 ■ 2^® - 61^ - 6 iGi(g)^ + 
o 1 y 413 

28iGi{g) - Hi - -Gi{g)^ + 8 Gi( 5 )" - y Gi( 5 )= - 2Gi{g) - —)g + 9 • 2®i" + 9 • 2HGi{g)^ - 27 • 2HGi{g) 

+ 1512®“°’^®'+! - 2712^®“°’^®^+! 
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-27 • + 81 • + 9 ■ 2®"^Gi(5)'‘ - 9 • 2<>~^Gi{gf - 9 • 2^~^Gi{gf + 15Gi(5)^2®“'^‘<®> - 27Gi(g)^2^®”‘^‘(®> 

+9 • 2®”^Gi(p) + 45Gi(g)2®”°*^®> - 81Gi(p)2^®"®*^®^ + 15 • - 27 • + 579 • 2®”^ 

+27 ■ 2^®+! - 5i' - 5iGi(g)2 + lbiGi{g) - 5i - ^Gi(g)" + ^Giigf + \Gi{gf - ^Giig) - 294]. 

(39) 


In this way, we have obtained an explicit formula for 
the ATT to an arbitrary node in Dg , as given by Eq. (l3^ . 
For the particular case z = 0 corresponding to the trap¬ 
ping with the trap placed at the central node, Eq. (15^ 
reduces to the expression in Eq. (I26E For another limit¬ 
ing case i = Mg that the trap is located at a peripheral 
node, we use Tp{g) to denote the ATT. For this case, 
plugging i = Mg = {g'^ — g + 2)/2 into Eq. (1391) yields 

Tpig) = r. ^ / o -^ fsi • 229 - 2^2 _ 189 . 2^3-2 

9-29-1-35-3 V y y 

189 • 2^9-1 - 27 ■ 2V + 81 • 23-^g'^ - 147 ■ 29-^5 
+123 • 29 + 3 / - y _ 58^2 _ 

(40) 

For very large systems {g —> 00 ), we have 

Tp{9) ^ Q 3 - - Ng{\nNgf , (41) 

which is also different from the Ng In Ng scaling of ATT 
for trapping in compact dendrimers with a peripheral 
node being the trap^i^. 

For general case of 0 < z < Mg, the expression for ATT 
provided in Eq. (15^ is rather lengthy and is hard to ob¬ 
serve the the dependence relation of ATT on the system 
size Ng. In order to show how the ATT scales with A’g, In 
Fig. [21 we report the results, given in Eq. (IMl) . for trap¬ 
ping in His with the trap positioned at nodes having dif¬ 
ferent distances towards to central node. Figure | 2 ] shows 
that the ATT depends on the level of trap node, which 
increases with the distance between the central node and 
the trap, i.e., Tr^{g) {i + l)Ng\nNg. Note that for 
trapping in compact dendrimers^, the ATT also grows 
with the level of the trap but displays a different behavior 
from that corresponding to extended dendrimers. 


C. Trapping with the trap uniformly distributed 

We are now in position to study trapping in ex¬ 
tended dendrimer Dg with the trap uniformly distributed 
throughout all nodes. In this case, the ATT, denoted by 
{T)g is defined as the average of MFPTs over all pairs of 
nodes in Dg : 




{T)s = 


- 1 ) 




(42) 


j=i 



FIG. 2. Average trapping time for the trap node at different 
levels in D 15 . 


For convenience, let Ttot{g) denote the summation term 
on the rhs of Eq. dUl), that is. 


Then, 


Ttotig) — Tij{g). 

i=i i=i 


{T),= 


Ttot {g) 


' NgiNg - 1 ) ■ 


(43) 


(44) 


By definition, {T)g involves a double average: The first 
one is the average of trapping times running over all start¬ 
ing nodes to a given trap node, the second one is the 
average of the first one with the trap having a uniform 
distribution. In what follows, we will apply two disparate 
approaches to analytically determine Ttotig) and {T)g, 
with the results obtained by both methods being consis¬ 
tent with each other. 

We first determine Ttot( 5 ) and {T)g by using the inter¬ 
mediary results obtained in Section IlIIBI For Ttot{g), it 
satisfies 


IVlg 

Ttotig) = E ^*( 5 ) [TrMiNg - 1)] . (45) 

i=0 

Inserting Eq. (|45ll into Eq. ([44|l and using Eq. ([39|, we 
obtain 

Mg 

j:N.ig)[TrMiN,-l)] 

{T)g = 


NANg - 1) 











1 

81.229-2 - 27 • 23g - 45 ■ 29-1 9^2 15^ g ' 

(^729 • - 1701 • - 1701 • 2®®"^ - 81 ■ 4®/ 

-1323 • 2®“^/ + 567 • 2^®”^^^ + 2835 • 2^®”^g 
-189 • 4® + 27 ■ 2®"^/ - 45 • 2®+^g^ 

-540 • 2®”^3^ - 207 • 2®g + 1233 • 2® + 39/ 

+129/ - 3/ - 7055 - 612^ . 

(46) 


Besides the above method, we can also compute Ttot{g) 
and {T)g by using the connection between effective resis¬ 
tance and MFPTs for random walks^i. For this purpose, 
we view Dg as an electrical networb^^ by considering ev¬ 
ery edge in Dg to be a unit resistor—. Let Rij{g) denote 
the effective resistance between a pair of nodes i and j 
in the electrical network associated with Dg. Then, the 
following relation holds 

"^ijig) + T^jiig) ~ 2^9 ^ijig) ■ ( 47 ) 

Using this link governing MFPTs and effective resistance, 
Eq. (H51) can be rewritten as 

JVb N, 

rtot(/=E,^^%(5). (48) 

i=l j=l 

Since extended dendrimers have a treelike structure, the 
effective resistance Rij{g) is exactly the usual shortest- 
path distance dij{g) between nodes i and j in Dg. Then, 
making use of Eqs. dm) and m, we have 

Ttot (g) = Eg dij {g) = Eg Plot (g) (49) 

1=1 


and 


/rp\ _ 7tot(g) 

^ “ Ng{Ng - 1) 


EgPtot jg) 
NgiNg - 1 ) 


- Eg ■ dg . 


(50) 


Inserting Eq. (ITSl) into Eq. (1501) . we recover Eq. (HOI) . 

When the system size is large enough, the ATT {T)g 
in Eq. (|40)) can be approximated as 


{T)g^ 


729 ■ 

81 • 229-2 


2®/~ TV/In TV/^ (51) 


which differs greatly from the leading asymptotical scal¬ 
ing for ATT in compact dendrimers with a trap dis¬ 
tributed uniformly among all nodes^. 


D. Result comparison and analysis 

From the above obtained results, provided by Eq. (1001) . 
Eq. (HIT) , and Fig. [OJ it can be seen that the ATT for 


trapping in extended dendrimers exhibits rich scalings 
when the trap is placed at different positions. Figure [0] 
and Eq. (l39l) show that the trap’s location has a sig¬ 
nificant effect on the efficiency of trapping in extended 
dendrimers measured by ATT; The trapping efficiency 
decreases when the distance from the trap to the central 
node increases. For the case that central node is the trap, 
the trapping process is the most efficient, its efficiency 
scales with system size Ng as Ng In Ng ; while for the case 
that a peripheral node is the trap, the trapping process is 
the least efficient, the efficiency varies with network size 
Ng as NgiluNg^. 

The different behaviors for (g) provided in Eq. (l39l) 
lie in the particular structure of the extended dendrimers. 
Figure [1] shows that the extended dendrimer under con¬ 
sideration consists of 3 branches (regions), each being a 
subtree with a node at level 1 as its root. For z = 0 cor¬ 
responding to the case that the central node is the trap, 
the walker, regardless of its starting position, will visit at 
most one branch before being absorbed by the trap. For 
T > 0 corresponding the case when the trap is placed on 
a node with distance i to the center, the walker, starting 
from a large group of nodes, must first visit the central 
node and then proceeds from the center along the path 

■Co — 'Cl — -- Ti-i — ri until it is trapped. Equation (1571) 

shows that the MFPT Tr^_^ri from a node at level z — 1 
to its direct neighbor at level z grows with z. Therefore, 
the scaling of Tr^{g) grows with increasing i. In particu¬ 
lar, when the trap is moved away from the central node 
to a peripheral node, the ATT changes from TVln(TV) to 
TVln(Af)2. 

On the other hand, when the trap is uniformly dis¬ 
tributed over Dg, the dominant scaling of ATT {T)g 
scales with the network size Ng as TVg(lnTVg)^, which is 
the same as that oi Tp{g) for trapping in Dg with a trap 
fixed at a peripheral node. Thus, the Ng In Ng behavior 
of the ATT to the central node is not representative of 
extended dendrimers, while the scaling of ATT to a pe¬ 
ripheral node is a representative property for trapping in 
extended dendrimers. The equivalence of the dominant 
scalings of ATTs for Tp{g) and {T)g can be easily under¬ 
stood. From Eqs. © and (nni), it is evident that in Dg 
the fraction of peripheral nodes is about 4. Furthermore, 
the ATT is the highest for all trapping problems in Dg 
with a single immobile trap. Therefore, the Tp[g) alone 
determines the leading behavior of {T)g. 

The aforementioned results also indicate that for the 
three cases of trapping problems under consideration, 
the trapping processes are less efficient in extended den¬ 
drimers than in compact dendrimers. Actually, trapping 
in extended dendrimers also displays different phenom¬ 
ena from other complex systems. For example, for trap¬ 
ping in Vicsek fractals^^*^! as a model of regular hyper- 
branched polymers22ri^, the efficiency is identical, de¬ 
spite of the location of the trap. Furthermore, trapping 
process in extended dendrimers is more efficient than in 
Vicsek fractals, implying that extended dendrimers have 
a desirable structure favorable to diffusion, as opposed 
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to Vicsek fractals. In addition to compact dendriiners 
and Vicsek fractals, trapping has also been studied in 
various other trees, such as T—shape treelike fractals 
and their extensions^^^— , fractal scale-free trees^, small- 
world uniform recursive trees^^— , non-fractal scale-free 
trees^. Previously reported results show that trapping 
behaviors in the first two network families are similar 
to those of Vicsek fractals, and that trapping behav¬ 
iors in the last two network families resemble those of 
compact dendrimers. Thus far, the dependence relation 
Ng{\nNg)^ of Tp{g) and {T)g on system size Ng has not 
been observed in other systems, implying that the struc¬ 
ture of extended dendrimers is unique. 

Before closing this section, we should stress that al¬ 
though we only focus on trapping in a particular ex¬ 
tended dendrimer, where any branching site has three 
nearest neighbors, for trapping in other extended den¬ 
drimers with the number of nearest neighbors of each 
branching site greater than 3, the qualitative scalings are 
similar to those found for the specific extended dendrimer 
under consideration. 


IV. CONCLUSIONS 

We have performed an extensive study of the trapping 
problem on extended dendrimers with varying lengths of 
the molecular branches. We first addressed a particu¬ 
lar case with the trap hxed at the central node; we then 
considered the case that the trap is positioned at an ar¬ 
bitrary node; and we finally studied the case with the 
trap uniformly distributed over all nodes in the systems. 
For these three cases of trapping problems, we studied 
the ATT as an indicator of the trapping efficiency, and 
obtained closed-form formulas for the ATTs, as well as 
their leading terms. We found that the location of trap 
has an essential influence on the trapping efficiency, with 
the ATT increasing with the shortest distance between 
the central node and the trap. We showed that the be¬ 
havior of trapping in extended dendrimers is unusual, 
since the leading scalings of the ATTs are unique, which 
are not observed in compact dendrimers and even other 
complex networked systems. 
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